function ng = uij (i,j)
  global Nx
  global Ny 
  ng = j + (i - 1)*Ny;
end

